A drop of hyperfine field at Sn in Fe/Cr/Sn/Cr multilayers 
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The magnetism of the Fe/Cr/Sn/Cr multilayers was studied by first-principles functional the- 
ory. The calculations by LAPW and SKKR methods showed that two solutions exist in the 
Feg/Cr-n/ Sn/Cr2 system. One of them is originated from the antiferromagnetic order in the 
bulk Cr, and the other is connected with incommensurate spin density wave (ISDW) in Cr. A char- 
acteristic feature of this system is realization of ISDW in the Cr film thinner than the half length 
period of a common ISDW in Cr. In the Cr layers with width lower than a quarter of the wave, 
the ISDW cannot be the cause of the sharp drop of the HFF at Sn. The calculations of the system 
Fes/Crg/Sn/Crg with and without allowance for the roughness at the Fe/Cr interface showed that 
the roughness leads to a significant decrease in the HFF at Sn. 

PACS numbers: 73.20.Dx, 75.70.Cn, 71.15.Ap, 71.15.Cr 



I. INTRODUCTION 

Since the discovery of the antiferromagnetic coupling 
between Fe layers p| and giant magnetoresistance in 
Fe/Cr multilayers a lot of studies have paid a great atten- 
tion to the Cr magnetic properties in these systems. Up 
to now, however, the discussion on the magnetic order 
in thin Cr films continues, there is no clarity in under- 
standing of the interconnection between the Cr magnetic 
state and the magnetic coupling of Fe layers. Recently, 
the authors of 3j succeeded in obtaining the Fe/Cr mul- 
tilayer systems with a Sn monolayer inserted. Multilayer 
systems Fe^/Cr^/ 'Sn/Crd 3 with different widths d\, di 
and c?3 were investigated by various physical methods in- 
cluding the Mossbauer spectroscopy. One of the most 
important results was the established dependence of the 
hyperfine magnetic field (HFF) at Sn nuclei on the width 
of the Cr layer. This dependence is characterized by a 
sharp drop of the HFF magnitude at a Cr width less than 
3 nm (d-2 + d% < 3 nm) 4] . Mossbauer experiments can- 
not give a direct answer as to the magnetization of Cr 
atoms, but if the HFF at Sn is assumed to be propor- 
tional to the local magnetic moments of the nearest Cr 
atoms, one should inevitably conclude that at a width of 
the Cr layer less than 3 nm the Cr magnetic moments in 
these systems are close to zero. 

As soon as the experimental data had been available, 
theoretical first-principles calculations of the electron 
structure and magnetic characteristics of Fe/Cr/Sn/Cr 
were conducted |3- 0- ■ The calculations were per- 
formed by both linear methods 0, IE Q and the Green 
function formalism 6] with different approximations of 
the exchange-correlation potential. As a whole, the re- 
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suits obtained by different methods and authors agree 
well with each other and do confirm the possibility of 
estimation of the Cr magnetization from Mossbauer ex- 
periments at Sn 6] . Some differences between the results 
have, though, led to difficulties in explanation of the pe- 
culiarities of the HFF variation and elicited a number of 
unsolved questions. 

In the paper presented, we concentrated our attention 
on the discussion of two hypotheses advanced earlier in 
0, Q for explanation of the Sn HFF behavior. The first 
one is connected with the existence of the high-spin and 
low-spin states and a transition between them with a 
change in the Cr width. Such solutions were obtained for 
the Feg/Cri4/Sn/Cr2 system in 0. The second hypoth- 
esis is connected with the effect of the Fe/Cr boundary 
roughness on the Cr magnetization. 

We have performed additional calculations of the peri- 
odical systems Fe§/Cru/Sn/Cr2 and Fe-ijCr^j SnjCr% 
(the latter being calculated with and without account 
of the roughness of the Fe/Cr interface). The calcu- 
lations were conducted by the following methods: the 
full-potential linearized augmented plane wave method 
(FP LAPW) realized in the program package WIEN2k 
8] , and a screened Korringa-Kohn-Rostoker method de- 
veloped in Julich [|| in the atomic sphere approximation 
(ASA SKKR) and in a full-potential scheme (FP SKKR). 
The calculations employ the exchange-correlation poten- 
tial in the local density approximation (LDA) in the for- 
malism of Refs. ^3 an( i E3- A s shown in Ref. if 
scaled to the moment, the LDA and the generalized gra- 
dient approximation (GGA) yield the same behavior. 

For our calculations we have used the following mod- 
els. The multilayer systems were presented as alternating 
layers in the (001) plane of a bcc lattice with two types 
of roughness (Fig.^ A and B) and without it. The inter- 
layer distances Fe-Fe, Fe-Cr, Cr-Cr and Cr-Sn were taken 
equal to 0.144 nm, which corresponds to the half lattice 
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FIG. 1: Two types of roughness at the Fe/Cr boundaries 
(solid lines). Dashed lines show the unit cells. 



parameter of the bulk Cr (0.288 nm). Notice that the 
lattice relaxation and small variations of the interlayer 
distances do not change the qualitative behavior of the 
results of Ref . . 



Results 

As stated earlier, the calculation of the 
Feg/Crxi/ SnjCr-i system by the SKKR method 
gave two solutions @ that we called "high-spin" (HS) 
and "low-spin" (LS) in accordance with the magnetic 
moment magnitudes (the difference in energy between 
the two solutions did not exceed 1 mRy per unit cell). In 
connection with these solutions it is appropriate to raise 
here some questions. First: why the LS solution was 
not found in Ref. [6j by the LAPW method? Second: 
do there exist two solutions in other systems and do 
they give a possibility to describe the peculiarities of 
the Sn HFF behavior as a transition from the HS to 
the LS state with a decrease of the Cr width? And the 
third question concerns the nature of the LS state (the 
character and peculiarities of the behavior of the Cr 
magnetic moments in the HS solution were considered 
in detail and explained in 



Two solutions were also found in an isolated 23-layer 
Cr film by other authors 13]. One of them was classified 
as a spin density wave and was connected with incom- 
mensurate spin density wave (ISDW) of bulk Cr. The 
ISDW has been studied from sixties, and now reliable ex- 
perimental and theoretical results are available for both 
bulk Cr, and some multilayer films with Cr, the majority 
of which are reflected in reviews 0, 0] in a rather full 
measure. We did not meet in literature with ISDW in 
Cr films thinner than 20-22 layers. This seems natural 
as the ISDW half period exceeds 20 layers. A hypothesis 
has, nevertheless, arisen that the nature of the low-spin 
state obtained in Q is connected with ISDW. For the 
existence of the ISDW solution, the Fermi-surface topol- 
ogy and the accuracy of its description are important. 
Assuming that in the FP LAPW calculation of the work 
the accuracy of the Fermi-surface description was un- 
sufficient, we have increased the number of basis wave 
functions limiting it by the parameter R.MTK max = 7.5 
(in it was RMTK max — 7.0). The expansion of the 
wave functions over the angular momentum was done 
up to l m ax — 10, and the Fourier expansion of the elec- 
tron charge density was conducted up to G max = 20. 
We should note that in the modification of the method 
used here, namely, the augmented plane wave plus local 
orbitals (APW+l.o.), the number of plane waves is usu- 
ally limited by R,MTK max — 7.0, as the convergence is 
achieved at lower number of basis functions than in the 
common realization of the LAPW method 0] . We have 
increased also the number of k-points in the irreducible 
part of the Brillouin zone up to 42, which may be crucial 
for the description of the Fermi surface topology. The 
calculation was conducted in a scalar-relativistic approx- 
imation for the valence electrons and it was fully rela- 
tivistic for the core electrons with an energy boundary 
at -7 Ry separating these two groups. 

The results of calculation of the Fe^/Cru/ Sn/Cr2 
system are given in Table [I] and Fig. [21 The calcula- 
tions were conducted with a charge convergence criterion 
10~ 4 -7- 10~ 5 (this is a squared charge-density difference 
between the last two iterations integrated over the space) . 
As seen in Fig. [2 in the FP LAPW method at the chosen 
parameters we have succeeded in obtaining two solutions, 
one of them being the LS solution. The HS solution is 
practically the same as in |(J by the FP LAPW method. 
The difference in energy between the HS and LS solu- 
tions is 1.4 mRy per cell. We should mention a very slow 
convergence of the LS solution, which prevented us from 
obtaining the accuracy over spin density of 10 -5 as for the 
HS solution (the accuracy of 10 -4 for LS was obtained). 
The LS solution obtained here has a peculiarity absent in 
the solution obtained by the SKKR method in Ref. [||, 
it is a phase slip between the 13-th and 14-th Cr atoms 
(if counted from the Fe layer, see Table [1}. The "phase 
slip" usually denotes a violation of the alternation of the 
magnetic moment directions characteristic for the anti- 
ferromagnetic ordering, and reveals itself in the fact that 
two neighboring Cr layers have equally directed magnetic 



TABLE I: The 

Fea/Cr 14 /Sn/Cr 2 . 



magnetic moments (ms) in 
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Fe 
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Fe 
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Fe 
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2.346 
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2.386 


Fe 
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17 
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Cr 


16 
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Sn 
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Cr 
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-0.007 


Cr 
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-0.014 


Cr 


12 
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Cr 


11 


-0.507 
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Cr 


10 
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0.157 


Cr 
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-0.531 
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Cr 
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0.284 


Cr 
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-0.343 


Cr 
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0.386 


Cr 
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-0.557 


-0.438 
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FIG. 2: The high- and low-spin solutions found by the FP 
LAPW method for the Fev/Cru/ 'SnjCri system 



moments. It is a common knowledge that the phase slip 
is one of the main features of the ISDW state 0]. To 
obtain an additional argument for identification of the 
LS solution as an ISDW state, we have recalculated the 
Feg/Cri4/Sn/Cr2 system by the SKKR method with an 
accuracy higher than in 

For the SKKR calculation a screening potential with a 



barrier height of 4 Ry is used and the structure constants 
include the coupling to the six neighboring spheres. A 
cutoff of the angular momentum at l m ax = 3 is used for 
the Green function, that implies a cuttof of the charge- 
density (in the case of ASA and FP) and potential (in 
the case of FP) components at 2l max = 6. The shape of 
the Voronoi cells for the FP calculations is described by 
the proper shape functions expanded up to 4i maa; = 12. 

The energy integration of the Green function for 
the determination of the electron density in the self- 
consistency procedure is performed by means of the 
Gaussian quadrature with 33 points on a semicircle (with 
5 Matsubara poles). For the k integration we use four 
kinds of meshes with 182, 45, 28, and 15 k points in the 
two-dimensional irreducible part of the Brillouin zone. 

The accuracy of the convergence over the potential was 
1CP 6 — 10~ 7 Ry, which is an order higher than in |6|. The 
increase of the number of points in these four kinds of 
&;-mesh up to 9800, 3150, 975 and 342, respectively, as 
well as the use of the double-dense energy mesh do not 
practically change the results. The calculation results are 
given in Fig. Eland 0] 
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FIG. 3: The ISDW solutions found by the ASA and FP 
SKKR methods 



From Fig. one may see that SKKR solutions also 
have a phase slip characteristic for the ISDW (between 
the 10-th and 11-th sites in ASA SKKR and between the 
11-th and 12-th sites in FP SKKR). Note that here we 
did not succeed in obtaining the second, HS, solution by 
the SKKR method with an accuracy better than in . 

Fig. 0] shows the calculation results at two tempera- 
tures (400 and 800 K). We should note that to describe 
the magnetic state variations with temperature (includ- 
ing the phase transitions), the Heisenberg excitations 
should be allowed for. Our calculations take into ac- 
count, with the help of the temperature Fermi distri- 
bution function, only Stoner-type excitations, which be- 
come effective at temperatures as high as 5000 — 10000 
K |l7j]. This allows us to suppose that, at 400 and 800 K 
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FIG. 4: The ISDW solutions found by the ASA SKKR 
method at temperatures 400 and 800 K 



used, the Stoner excitations weakly affect the magnetic 
moment magnitude, and to consider the temperature as 
a mean to change the details of the Fermi surface that 
determine the position of the phase slip. Indeed, from 
simple considerations, we may suppose that the ISDW 
period is determined by the position of the chemical po- 
tential which depends on the temperature (see, for ex- 
ample, [3) 

fi = €p — n 2 ksTn '(ep) /6n(e_p) + o{ksT / cf) 

where ep is the Fermi energy, ks is the Boltzmann con- 
stant, n{eF ) and n'(ep) are the electron density of states 
at the Fermi level and its derivative (the value of n'(ef) 
is negative in our system). With increasing temperature, 
the chemical potential grows, decreasing thus the dis- 
tance in the k-space between the electron and hole sur- 
faces which determines the ISDW period. That is why 
the phase slip moves from the 10-11-th sites to the 11-12 
sites with the temperature change from 400 to 800 K (see 
Fig. EJ. 

The above results leave no doubt that the nature of 
the LS solution in the Feg/Cru/ 'SnjCr-i system is con- 
nected with the ISDW state in Cr. The realization of this 
state in our system has, however, some features absent in 
the ISDW in other systems |14.llq| . The Sn layer present 
in the system gives the possibility to form a state with a 
quarter of the ISDW, which was not observed in theoret- 
ical calculations earlier. Furthemore, the role of the Sn 
layer in the ISDW formation is here different from that 
found in [19j in the Cr/Sn multilayers, where Sn layers 
are pinning centers for the spin-density wave maxima. 

Proceeding from the nature of the LS state one may 
assert that its realization in the multilayer systems with 
a Cr layer width less than a quarter of the ISDW (11-12 
Cr layers) is impossible, so the hypothesis about the exis- 
tence of the two solutions (HS and ISDW) in thinner Cr 
layers and the transition between them with a decrease 



TABLE II: The HFF at Sn (T) and magnetic moments (fi B ) 
in Fe 3 /Cr s /Sn/Cr$ 



Atom 


layer 


ideal 


A roughness (see Fig. 


B roughness 


Sn HFF 




-16.1 


-5.0 


-4.7 / -5.0" 


Sn 


11 


0.01 


0.004 


0.004 / 0.004° 


Cr 


1 


0.73 


0.23 


0.23 


Cr 


2 


-0.48 


-0.15 / -0.15" 


-0.15 / -0.15 a 


Or 


Q 
O 


U.OU 


U.lo 


U.lo 


Cr 


4 


-0.57 


-0.18 / -0.18" 


-0.16 / -0.18 a 


Cr 


5 


0.60 


0.19 


0.17 


Cr 


6 


-0.55 


-0.18 / -0.24" 


-0.15 / -0.22 a 


Cr 


7 


0.54 


0.23 


0.19 


Cr / Fe 


8 


-0.60 


-0.53 / 1.92 a 


-0.50 / 1.89" 


Fe 


9 


1.97 


2.25 


2.25 


Fe 


10 


2.46 


2.43 


2.37 / 2.47 a 



.../... denotes values for two nonequivalent atoms in the layer 



of the Cr width is untenable. This is confirmed also by 
calculations of the systems with thin Cr layers |6j , where 
no low-spin solution was found. 

So, the second explanation of the experimental depen- 
dence of the HFF at Sn connected with the effect of the 
Fe/Cr boundary roughness becomes more preferable. To 
check this hypothesis we have considered two models of 
the roughness of the Fe/Cr boundary shown in Fig. ^ 
The calculations for the film Fe^/Crs/ Sn/Cr$ were con- 
ducted by the FP LAPW method, the results are given 
in Table ITT1 One may sec from Table [H] that the rough- 
ness in both model cases leads to a decrease in magnetic 
moments of the Cr atoms nearest to the Sn layer. Tis 
results in a significant decrease in HFF at Sn from 16.1 
T (for the system with ideal interfaces) to 4.7 - 5.0 T 
(for the systems with roughness). This does not explain 
in full measure the experimental sharp drop of the HFF. 
To obtain it, one should perform the calculation of the 
systems with higher width and different roughnesses. 



Summary 

The first-principles FP LAPW calculation of the mul- 
tilayer model system Feg/Cru/ SnjCr-i gives two solu- 
tions. One of them, with a lower magnitude of the Cr 
magnetic moments and a phase slip, was identified as an 
ISDW. This identification was confirmed by the SKKR 
calculations. A characteristic feature of the ISDW state 
in our case is the fact that the Sn layer promotes the 
realization of the ISDW state in the films with Cr width 
less than the half period of ISDW. 

In connection with the above statements, one may as- 
sert that the band calculations of the periodical multi- 
layer systems without the interface roughness cannot ex- 
plain the experimentally observed sharp drop of the HFF 
at Sn with a decrease of the Cr layer width. 

The presence of roughness in the Fe/Cr interfaces leads 
to an HFF drop in thin films. This is confirmed by the 
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first-principles FP LAPW calculations conducted for the 
model Fes/Cr s / Sn/Cr% with ideal Fe/Cr interfaces and 
with two types of roughness. The calculations showed 
that there is a threefold drop in HFF in the systems with 
roughness as compared to the system with ideal Fe/Cr 
boundaries. However for explanation of the experimental 
dependence it is necessary to calculate a number of sys- 
tems with more than 20-23 Cr layers and with allowance 
for the different roughness. 
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